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MEASURING POSITION AND O RIENTATION USING MAGNETIC FIELDS 
Background of the Invention 
The present invention relates to a system and method 
for determining the position and orientation of a remote 
object relative to a reference coordinate frame using 
magnetic fields. More specifically, the position and 
orientation of a medical device, such as a catheter, within 
a patient is determined. 

Various magnetic tracking systems use circular loop 
antennas for magnetic field generation and detection. One 
or more coils generate distinguishable magnetic fields and 
one or more coils detect the generated magnetic fields. 
Knowing the position, orientation and geometry of the 
generating coils and the signal sensed in response to the 
magnetic field generation, the position and orientation of 
the sensing coils can be determined with respect to the 
field generators. The magnetic field generators usually 
establish the coordinate reference frame from which 
measured position and orientation is referenced. Depending 
on the number of coils used for generation and sensing, 
measurement of values for up to 6 degrees of freedom (three 
position values and three orientation values) of an object 
can be obtained. The equations normally used for 
calculating position and/or orientation using a generated 
magnetic field are almost all based on the assumption that 
the distance between the generating and detecting coils (r) 
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is larger than the radius (a) of the generating coil. This 
equation is known as the infinitesimal dipole model. This 
assumption is explicitly used in the following U.S. 
patents, hereby incorporated by reference. 

Hansen (4,622,644) discloses a 5 degree-of-freedom 
(5DOF) position and orientation (P&O) measurement system 
that uses 2 or 3 sets of orthogonal Hall effect sensors and 
a bar magnet for generating a magnetic field. It is noted 
that using 2 sets of 3-axis orthogonal Hall effect sensors 
restricts the region of operation due to ambiguities. a 
noise error vector is noted that contains the estimated 
residual noise measured at the sensors. His method of 
solution starts by determining an approximate starting 
point for a calculation process that solves the system of 
15 9 non-linear equations in 5 unknowns (x position, y 
position, z position, bar magnet azimuth, bar magnet 
elevation). The non-linear equations are based on the B- 
field approximation. Likewise, Hansen (4642786) uses the 
same P&o solution approach in a variety of tracking system 
20 configurations (both 5D0F and 6D0F) using the 
retransmission of magnetic fields. 

Rorden et al. (4710708) discloses locating methods 
using various generating and detecting coils for 
underground surveying. These include 2, 3, or more, 3-axis 
orthogonal detecting (sensing) coil sets with a single 
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transmitting coil; a 2-axis orthogonal generating coil set 
with a 3-axis orthogonal detecting set; a plurality of 
single axis collinear generating coils with a three axis 
orthogonal detecting set, among others. It is noted that 
5 at least 6 coil elements are required to avoid ambiguities 
or additional information must be utilized. Also noted is 
that using 2 sets of 3 axis orthogonal detecting sets 
restricts the region of operation due to ambiguities. The 
method of solution utilizes a process that solves the 
10 system of 6 non-linear equations in 6 unknowns (x position, 

y position, and z position, x direction cosine, y direction 
cosine, and z direction cosine). The direction cosines are 
composed of trigonometric combinations of azimuth and 
elevation such that this set of equations can be reduced to 
15 5 unknowns. The non-linear equations are based on the B- 
field approximation. 

Bladen et al. (WO 94/04938) discloses two system 
configurations: 3 3-axis orthogonal generating coil sets 
with a single axis detecting coil (5DOF) and one 3-axis 
20 orthogonal generating coil set with a 3-axis orthogonal 
detecting coil set (6DOF) . Bladen' s Appendix A describes 
the approximate B-field equations utilized in their 
process. Like Hansen, the 5DOF method of solution starts 
by determining an approximate starting point for the 
25 process and then solves a system of 9 non-linear equations 



- 3 - 



WO 99/32033 



PCT/US98/13431 



m 5 unknowns (x position, y position, z position, azimuth 
and elevation) . 

Dumoulin et al. (5377678) starts with equations that 
are remarkably similar to Rorden' s equations (again, B- 
field approximations). He never mentions the magnetic 
moment (m) from Rordens equations and has a 5DOF system 
with 5 receiver coils and a single axis generating coil. 
He solves 5 non-linear equations in 5 unknowns (x position, 
y position, z position, azimuth and elevation) and notes 
that there must be at least 5 coil elements. To start the 
process he places the generator coil at a known position 
and orientation. As an aside, he first states that two 
transmitter coils are needed to determine both position and 
orientation (column 2, lines 66 on) and then states he 
needs one transmit coil to determine position and 
orientation (column 4, lines 41 on). 

While it may not be obvious, a number of other patents 
use the infinitesimal dipole approximation, including Jones 
(4737794, 5307072), Raab (4054881, 4314251, 4346384), Crow 
et al. (4328548), Egli et al. (4287809, 4394831), Constant 
(4396885), Blood (4613866) and Kuipers (4742356). 

Published PCT application WO 97/32179 to Acker 
discloses a location system wherein at least one field 
generating coil or device creates different electromagnetic 
fields during different time intervals. This system 
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specifically covers pulsed-DC field generation and 
excitation control. It also describes a non-linear 
equation solution to the tracking problem, which has 
already been covered by others. 

A few patents note needed corrections when using the 
infinitesimal dipole model to make it more accurate when 
working close to the generating coil, Jones (4737794) 
builds upon the dipole model and discloses a quasi-closed 
form (non-iterative) matrix solution of P&O using one 3- 
axis orthogonal generating coil set with a 3 -axis 
orthogonal detecting coil set (6DOF) . As noted above, 
Jones starts with the dipole approximation. He then 
introduces additional correction terms from an infinite 
series of Legendre polynomials with the unstated proviso 
that the distance between the generating and receiving 
coils (r) is greater that the transmitting coil radius (a) . 
These correction terms are calculated based on the previous 
P&O solution. He casually mentions that elliptic integrals 
can also be used for generating these correction terms to 
the dipole model. 

Jones (5307072) discloses methods for compensating for 
positional offsets of generating coil and sensing coil 
sets. These methods are deeply embedded into his P&O 
process disclosed in 4737794. As noted above, Jones starts 
with the dipole approximation. He then introduces 
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additional correction terms from an infinite series of 
Legendre polynomials, again with the unstated proviso that 
r>a that correct for coil size, geometry and coil 
placement. These correction terms are calculated based on 
the previous P&o solution but can be used directly. The 
use of elliptic integrals, hypergeometric functions, etc. 
is also mentioned as a means of generating these correction 
terms to the dipole model. The methods detailed in both 
Jones patents will yield incorrect results when r<a. 
Neither will work at all when r=0 (the sensor coil 
collocated within the 1 generating coil). 

Extending these concepts further, certain patents 
address generating coils of any shape. They also require 
the use of non-linear equation solving to provide part or 
15 all of their position and orientation solution. Ben-Haim 
et al. (WO 96/05768) discloses a position and orientation 
measurement system that uses any shape generating coils. He 
too uses Legendre polynomials to define the magnetic field 
structure, correct for r>a only. He claims as new to the 
art non-concentric sensing coils for his 6DOF approach but 
Jones (5307072) has already addressed this issue. 

Blood (5600330) discloses a system to determine P&O 
using 2 or 3 non-dipole generating antennas (e.g., 
rectangular coils) that are not necessarily orthogonal. It 
does use a 3 axis orthogonal detecting coil set and yields 
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a 6DOF solution. The distortion immunity and large volume 
of measurement space are obtained by constraining the 
sensing device to work within the perimeters of the 
generating coils. it is noted that a process that 
iteratively determines P&O could have been developed. It is 
also noted that the number of receiving antenna (i) times 
the number of transmitting antenna (j) must equal at least 
6 with the requirement that i and j must at least be two 
(for 6DOF). It is noted that the location, orientation and 
geometry of the coils must be known, implying that the B- 
field equations can be explicitly stated. A non-linear 
process solving 6 equations in 3 unknowns (x position, y 
position and z position) based on measured magnitudes and 
vector dot products is developed. An exact non-iterative 
process is then used to compute orientation. Bloods 
process eliminates the need to compute elliptic integrals 
or Legendre polynomials but only works well within the 
confines of the generating coils. 

Other patents that use non-dipole fields in very 
different ways include Acker et al. (5558091), Martinelli 
(5592939) and Voisin (5172056). These disclosures use 
multiple coils to form fields that are not dipole in nature 
but vary quasi-linearly in certain respects. 
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pb i ects Summary af the Invent inn 
Accordingly, it is a primary object of the present 
invention to provide a new and improved position and 
orientation ( P&0 ) system allowing determination of the 
5 position and orientation of a remote device. 

A more specific object of the present invention is to 
provide a new and improved medical system. 

A further object of the present invention is to 
provide a laser catheter with precise position 
10 determination. 

Yet another object of the present invention is to 
provide a catheter having position determination without 
requiring multiple bulky structures in the catheter. 

A further object of the present invention is to 
provide catheter position determination without placing 
constraints on the relative positions of the catheter and 
an external locating system. 

Yet another object of the present invention is to 
provide a catheter for percutaneous myocardial 
20 revascularization. 

The above and other features of the present invention 
which will be more readily understood when the following 
detailed description is considered in conjunction with the 
accompanying drawings are realized by a medical system for 
25 treating a patient. Pluralities of transducers are 
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operable to determine the location of a medical device in 
a patient, the transducers including one or more device 
transducers and one or more external transducers. A 
medical device is insertable into a patient, the medical 
5 device having the one or more device transducers. An 
external locator is operable outside and adjacent to the 
patient, the external locator having the one or more 
external transducers. The system sequentially energizes 
one or more of the device transducers and/or the external 

10 transducers. Sensed signals are generated from one or more 
of the transducers, each sensed signal from a given one of 
the transducers depending on the energizing of another of 
the transducers. A processor is responsive to the sensed 
signals, the processor determining at least two location 

15 parameters of the medical device in the patient using a 
plurality of equations, which equations do not require a 
minimum separation distance from the energized transducer 
to the transducer generating the sensor signal for 
validity. As used herein, location parameters shall include 

20 position and orientation information relative to the 
catheter (or other remote device) . 

The present invention may alternately be described as 
a method of determining the location of a remote device 
comprising the steps of: using a system having a plurality 

25 of transducers operable to determine the location of a 
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remote device, the transducers including one or more device 
transducers and one or more locator transducers; placing 
the remote device in an operation zone, -the remote device 
having the one or more device transducers; placing a 
locator away from, but within an interaction distance from, 
the remote device, the locator having the one or more 
locator transducers; sequentially energizing one or more of 
the device transducers and/or the locator transducers; 
generating sensed signals from one or more of the 
transducers, each sensed signal on a given one of the 
transducers depending on the energizing of another of the 
transducers; and supplying the sensed signals to a 
processor, the processor determining at least two location 
parameters of the remote device using a plurality of 
equations, which equations do not require a minimum 
separation distance from the energized transducer to the 
transducer generating the sensor signal for validity; and 
wherein the system is a medical system used to treat a 
patient, the remote device is a medical device and the 
operation zone is within a patient. 

The present invention may alternately be described as 
a method of determining the location of a remote device 
comprising the steps of: using a system having a plurality 
of transducers operable to determine the location of a 
remote device, the transducers including one or more device 
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transducers and one or more locator transducers; placing 
the remote device in an operation zone, the remote device 
having the one or more device transducers; placing a 
locator away from, but within an interaction distance from, 
the remote device, the locator having the one or more 
locator transducers; sequentially energizing one or more of 
the device transducers and/or the locator transducers; 
generating sensed signals from one or more of the 
transducers, each sensed signal on a given one of the 
transducers depending on the energizing of another of the 
transducers; and supplying the sensed signals to a 
processor, the processor determining at least two location 
parameters of the remote device using a plurality of 
equations, which equations do not require a minimum 
separation distance from the energized transducer to the 
transducer generating the sensor signal for validity; and 
wherein each of the transducers generates and/or senses an 
electromagnetic field. 

Brief Description of the Drawings 
The above and other features of the present invention 
will be more readily understood when the following detailed 
description is considered in conjunction with the 
accompanying drawings wherein like characters represent 
like parts throughout the several views and in which: 
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FIGS. 1-5 are simplified xyz schematic diagrams for 
explaining the derivation of various equations and 
configurations used with the present invention; 

FIG. 6 is a front view of a two coil structure used 
with some embodiments of the present invention; 

FIGS. 7A, 7B, and 8-10 are simplified xyz schematic 
diagrams for explaining the derivation of various equations 
and configurations used with the present invention; and 

FIGS. 11 and 12 are schematics of two different 
embodiment systems according to the present invention with 
5 degrees of freedom (DOF) measurement capability. 

FIGS. 13 and 14 are schematics of two different 
embodiment systems according to the present invention with 
6DOF measurement capability. 

FIG. 15 is a simplified top view of a table and gimbal 
arrangement used for calibrating the present system. 

FIGS. 16 and 17 are flowcharts on the use of the norm 
and residuals comprising the norm as a safety mechanism and 
a method of compensating for magnetic distortions. 

FIG. 18 is a flowchart on the use of a known starting 
point as a means of calibrating a coil in a catheter. 

FIGS. 19-23 are simplified top views of various 
disposable sensors used with the present invention. 
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Detailed Description 
A number of magnetic position and orientation 
determining system configurations using magnetic fields are 
presented. The methods determine the relative position and 
orientation of a sensing device with respect to a generator 
of magnetic fields which establishes the reference 
coordinate frame. The sensing devices can be one, two or 
three-axis magnetic sensing elements such as coils of wire, 
Hall-effect sensors, magnetometers and/or any sensors based 
on magnetoresistive, magnetoinductive, or galvanomagnetic 
technology. The field generating devices can be one, two 
or three-axis coils of wire or magnets. All configurations 
can be reversed using the principal of reciprocity; i.e., 
the sensing and generating elements can be interchanged. 
The systems can utilize AC waveforms or pulsed-DC waveforms 
or retransmitted waveforms. Time, frequency, phase or 
various combinations of each can be utilized to distinguish 
between the various excitations of the generating elements. 
The process can also be used to improve the solutions of 
other tracking systems that are based on approximate B- 
field formulations. 

All the configurations calculate the position and 
orientation solution based on the exact formulation of the 
near field (or quasi-static) magnetic field coupling as 
opposed to approximations used elsewhere. Of course, a P&O 
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process utilizing equations for the far field (radiation) 
magnetic field could also be solved. Generating coil 
position can be included in the equation by making 
substitution of the form: 
x' =x-xo 
y' =y-yo 
z' =z-zo 

where x' , y' and z> are the translated positions from the 
true x, y and z coordinate system with offsets xo/ yo and 
ZO. Generating coil orientation can be included in the 
equations by multiplying a 1x3 B-field vector with a 3x3- 
rotation matrix. 

The process utilizes the known position, orientation 
and geometry of the generating elements, their effective 
areas, driver gains, preamplifier gains, sensor effective 
areas and other system constants. The effective area for 
a coil is the area of a single coil turn times the number 
of turns. This calibration data is gathered during the 
manufacturing process and is applied in real time during 
execution of the P&o process. 
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Two examples of exact B-field equations are: 

Circular Coil of Wire 
Coil in the x-y plane (figure 1): 



XT 



B 



W y 



a 1 +x' +y* 



-K(a)+£(a) 



B. 



2n 



K(a)+£( a ) 



(2) 



(3) 



where: 



10 



*(«) = I 



"4 



(4) 



£(a)= /VO -a'an J 8 )/9 (5) 
p 

j 4 or 



(6) 



H Q is the permeability of free space; a is the coil radius; 
I is the current flowing through the coil and X represents 
the number of turns comprising the coil. Equations for 
coils in the x-z plane and y-z plane can be obtained from 
the above by the use of obvious substitutions or by using 
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a 3x3 rotation matrix. 

Rectangular Coil of Wire 
Coil in the x-y plane (figure 2): 



B - X ^ 7 (-b+2y)z 
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M 0 is the permeability of free space; a is the length of the 
loop in the x direction; b is the length of the loop in the 
y direction; I is the current flowing through the coil and 
A represents the number of turns comprising the coil. 
Equations for coils in the x-z plane and y-z plane can be 
obtained from the above bythe use of obvious substitutions 
or by using a 3x3 rotation matrix. 

Any means of generating a B-field that can be 
explicitly stated in mathematical terms can be utilized in 
this patent. it should also be noted that mathematical 
approximations to the B-field will work, but only with as 
much accuracy as the approximation allows, e.g., an 
approximation based on the infinitesimal dipole only works 
for r»a. Some of these approaches are already patented. 

Some words about the rectangular vs. circular coils, 
are in order. The evaluation of magnetic field generated 
by rectangular coils need not be expressed in terms of 
elliptic integrals. While the exact form of the B-field 
for rectangular coils appears foreboding, the elimination 
of the need for calculating the values of the elliptic 
integrals may in fact be a plus. (Elliptic integral 
evaluation is usually based on rational function 
approximations or other methods that take extra time to 
compute.) This is especially true if the Jacobian (or 
Hessian) needs to be calculated (see following section). 
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Exact calculation of the Jacobian and/or Hessian can be 
accomplished using symbolic software. The Jacobian and 
Hessian for the circular coils, however, require further 
evaluations of the elliptic integrals. 

While it is not my intent to teach others about 
solving non-linear systems of equations or least squares 
(optimization) problems, a few words are in order. There 
are two types of non-linear equation problems, the non- 
linear systems of equations problem and the non-linear 
least squares problem. The non-linear systems of equations 
problem is: given n functions f 1# f 2 ...,f ft of the n variables 
x 1f x r ,x n find values for x 1 ,x 2 ...,x n that solve the non-linear 
system of 

f,(x v x 2 ...,x n )=0, l<i<n 

The non-linear least squares problem is: given m functions 
f 1# f 2 ,... / f BI of the n variables x 1# x ? ...,x n , with m>n, find 
values for x 1 ,x 2 ...,x n that solve the non-linear least squares 
problem 
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Many methods exist to tackling these problems with the 
predominant methods requiring the evaluation of the 
Jacobian (a matrix of partial derivatives of the equations 
with respect to the unknowns), either explicitly or by 
finite differences and sometimes requiring the evaluation 
of the Hessian (a matrix of second partial derivatives of 
the equations with respect to the unknowns) . These methods 
are often referred to as Newton methods, gradient methods 
or steepest descent methods or variations on that theme. 
The discussion that follows refers to various publications 
listed and numbered per the list of documents incorporated 
near the end of this application. References l, 2 and 11 
are good introductions to these methods and provide basic 
processes. References 3 and 5 provide more advanced 
versions of these methods. Because these methods are 
iterative there is no guarantee that they will converge to 
the correct solution and convergence is only guaranteed 
when the starting point for the iteration is 'close" to the 
final solution. This implies that these methods all need 
a starting point for the process. This will be discussed 
in more detail under 'Starting Point for Process.' Most of 
these methods incorporate techniques to expand the region 
of convergence (so called global methods). Non global 
methods are the ones that have been noted in Hansen 
(4622644, 4642786), Rorden et al. (4710708), Dumoulin et 
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al. (5377678) and Blood (5600330). It should be noted that 
increasing the region of convergence is most important when 
first starting the process. It allows the process to 
converge from poor starting points. Once started, however, 
a limited region of convergence is ok since the next 
starting point is based on the last solution that is 
usually close to the present solution. 

Reference 11 also discusses the method of simulated 
annealing, a controlled random search method for least 
squares problems. Reference 8 is another statistically 
based search method. Reference 4 and 12, useful for 
solving non-linear system of equations, are based on 
successive substitutions instead of simultaneous solutions 
(as done in Newton methods) . They are derivative free and 
are related to Gaussian elimination in linear equations. 
Reference 7 is a method of formulating the solution of the 
non-linear system of equations into a system of ordinary 
differential equations. Standard differential equations 
processes are than used to solve the problem. Reference 9 
and 10 are bisection methods that attempt to find the roots 
of a system of non-linear equations by bisecting an n 
dimensional space in a systematic manner. 

Homotopy/continuation methods are discussed in reference 
13. These methods work by constructing a continuous 
deformation of one set of non-linear equations into the 
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10 



original system of non-linear equations. Reference 14 
covers genetic processes that try to evolve a best 
solution. 

All these methods (and others not mentioned) can be 
used to solve the P&o processes that are to be discussed. 
The various methods require different processing power 
depending on the size and complexity of the problem. In 
general, bisection, homotopy, genetic and statistically 
based methods take much longer to converge than the Newton 
based methods and are therefore not feasible at the present 
time for real time use. Bisection and homotopy methods, 
however, can usually find the global solution from any 
starting point. Newton methods can be used for real time 
computation and global Newton methods extend the range 
between the starting point of the process and the final 
solution. The methods based on Gaussian elimination are 
the fastest of the lot but require a starting point 
relatively close to the solution. The best strategy for 
starting the process at an arbitrary point might be to use 
one of the slower methods at start up and switch over to 
Newton or Gaussian methods once convergence has been 
achieved. 

A point worth mentioning is that the set of equations 
for the different tracker configurations is variable and 
25 that a particular set might be better based on the 
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placement or configuration of the coils. That is, if there 
are more equations than unknowns, one can formulate the 
problem as: 

1) Solve a system of non-linear equations (m=n) . 

2) Solve the non-linear least squares problem (m>n) . 
Solving a system of non-linear equations can imply using a 
subset of the equations. Using the least squares approach 
might mean using all the equations or a subset of the 
equations. Those familiar with the problem will obviously 
be able to put together other sets of equations not 
mentioned here. Using more than the minimum required 
number of equations using the least squares approach 
usually reduces the sensitivity of the P&O solution to the 
model parameters and increases the range of convergence. 

Up to this point all discussion of the exact 
formulation of the underlying equations has been directed 
towards the B-field generating elements. However, sensing 
devices, such as coils, are of finite size. The use of the 
exact equation formulation for the sensed voltages might be 
necessary depending on the applications accuracy 
requirements and the sensor coil size. As has been noted 
in most of the applicable patents the roles of the B-field 
generators and the B-field sensing elements can be 
reversed. A similar approach to correcting for the finite 
size of the sensing coils can be derived and utilized in 
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the P&O process. For most applications, especially those 
in catheters, this added complexity to the defining 
equations is neither required nor desired. 
Configurations 

All of these configurations can be written in a 
general mathematical format. The general formulation is 
(in terms of an orthogonal dot product): 

•A.fcoia r.coiP y.co*r-*i*v„ (10) 

Or, in terms of the problem formulation, 

[b.J, J*\ # K t os a ,T ' C0 * P • J. co * T •* 1 - ■ D (M) 

where B^, B y and B 2 are the x, y and z components of the p th 
generating coil' s B-field; i, j and k are unit vectors in 
the direction of the Cartesian coordinate axes; cos(a), 
cos(0) and cos(y) are direction cosines of sensor coil q; 
k q is the effective area of sensing coil q; k ? represents 
the number of turns comprising the p th generating coil and, 
for time varying currents, a term proportional to the 
negative time derivative of the B-field; and v ft m is the 
induced voltage at sensing coil q for the p th generating 
coil. The unknowns are the x, y and z position of the 
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sensing coil and the 3q direction cosines. For 1 sensing 
coil q=l . 

It is also possible to rewrite the above in terms of 
the sensing coils azimuth (A) and elevation (E) . cos(a), 
cos(0) and cos(y) are then replaced by cos (A) cos (E) , 
sin(A)cos(E) and -sin(E), respectively. The unknowns are 
now one less than before, i.e., x, y and z position of the 
sensing coil and sensing coil azimuth and elevation (hence 
5DOF) . 

4 Generating Coils, 1 Sensing Coil 

4 generating coils (p=4) can be placed at arbitrary 
positions and orientations as long as they are not 
collocated and oriented the same way. The coils should be 
arranged such that symmetries are minimized to prevent 
places having the same B-fields for different positions and 
orientations. The 4 generating coils may be coplanar or be 
two sets of orthogonal coils or anything else (figure 3). 
The 4 B-field equations can be solved as a non-linear 
systems of equations problem by adding the additional 
equation that the sum of the squares of the direction 
cosines must sum to 1, i.e., 

cos 7 a + cos a p + cos 2 y - 1 = o (12 ) 
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which is equivalent to 

(cosAcosE) 1 +(sin^cos£) 1 + (-sin£) 2 -1 = 0 (13) 

Now m=5 and n=5. 

5 Generating Coils, 1 Sensing Coil 

5 5 generating coils (p=5) can be placed at arbitrary 

positions and orientations as long as they are not 
collocated and oriented the same way. The coils should be 
arranged such that symmetries are minimized to prevent 
places having the same B-fields for different positions and 

10 orientations. The 5 generating coils may be coplanar or be 
two sets of collocated orthogonal coils with a fifth coil 
arbitrarily placed or anything else. The 5 B-field 
equations can be solved as a non-linear systems of 
equations problem (m=5, n=5) or solved as a least squares 

15 problem by adding the additional equation that the sum of 
the squares of the direction cosines must sum to 1, as 
above. Now ra=6 and n=5. 

6 or More Generating Coils, 1 Sensing Coil 
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6+ generating coils (p=6+) can be placed at arbitrary 
positions and orientations as long as they are not 
collocated and oriented the same way. The coils should be 
arranged such that symmetries are minimized to prevent 
places having the same B-f ields for different positions and 
orientations. The 6+ generating coils may be coplanar or 
be two sets of 3-axis orthogonal coils or 3 sets of 2-axis 
orthogonal coils or anything else. The 6 coil problem can 
be solved as a non-linear systems of equations using the 6 
B-field equations (ro=6, n=6) where the unknowns are x, y 
and z position and the 3 direction cosines. Or it can be 
solved as a least squares problem by adding the additional 
equation that the sum of the squares of the direction 
cosines must sum to 1, as above. Now m=7 and n=6. Adding 
additional coils adds additional equations (increasing m) 
for the least squares problem to fit to. While this can 
reduce the errors in the P&O solution and increase the 
region of convergence, it does so at the expense of 
processing time. 

A special configuration of generating coils should be 
noted. Hansen (4622644) and Rorden et al. (4710708) both 
note that there are ambiguities when 2 sets of 3 collocated 
orthogonal generating coils are used. This can be best 
understood by examining a specific configuration where the 
2 sets of triads are located symmetrically about the origin 
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on the y-axis at a distance y 0 as shown in figure 4 • The 
sensing coil is place exactly between the two triads (y=0) 
and is oriented such that cos(a)«=l, cos(0)=o and cosfY) 5 ^. 
Writing the equations one sees that there are only 4 unique 
equations, i.e., equation (14) and (15) are identical with 
equation (17) and (18). This results in a matrix that is 
not of full rank (there are only 4 independent equations 
from the set of six) . The mathematical methods used by 
Hansen and Rorden et al. could not handle this problem. 
However references 3 and 5, among others, can solve rank 
deficient systems of non-linear equations and least squares 
problems without encountering the problems noted by Hansen 
and Rorden et al. 



n 



(14) 



(13) 
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(17) 



Vgy ■ A 4 A, 



(it) 



(19) 



Also, as noted earlier in the 4 Generating Coils, 1 Sensing 
Coil configuration, it is possible to add the sensing coil 
normal equation to this configuration and obtain 5 
equations with the 5 unknowns of X, y, z, azimuth and 
elevation. 

2 Generating Coils, 2 Sensing Coils 



2 generating coils (p= s 2) can be placed at arbitrary 
positions and orientations as long as they are not 

10 collocated and oriented the same way. The coils should be 
arranged such that symmetries are minimized to prevent 
places having the same B-fields for different positions and 
orientations. The 2 generating coils may be coplanar or 
collinear or be a set of collocated orthogonal coils or 

15 anything else. 
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The 2 sensing coils (q= 2 ) can be orthogonal, coplanar 
or collinear or some other combination as long as they are 
not collocated and oriented the same way (figure 5). There 
are now 2 B-field equations per 2 sensing coils for a total 
of 4 equations. A 5DOF system is obtained if the 2 sensing 
coil normals are collinear. Coplanar coils might seem to 
yield a 5DOF system since the coils are both oriented the 
same way. However, 6DOF is possible because the 2 coil 
centers define a plane without axial symmetry. Rotation 
about the sensing coil plane' s x and z-axis is defined by 
the relative positions of the 2 coils. This and the coil 
normal are sufficient to determine the 3DOF required for 
orientation. For the case where they are not oriented in 
the same direction, e.g., orthogonal, a 6DOF system is also 
obtained. Expressed in terms of equation (11), there are 9 
unknowns. The unknowns are the x, y and z position of the 
sensing coils and the 6 direction cosines. 

Two or more additional equations are still needed to 
use the methods mentioned above. Two obvious equations 



are: 



2 2 5 

cos a, + cos * P J[ + cos *y - 1 = o (20) 

cos 2 aj + cos 2 |3 J + cos J y J -l= : 0 (21) 
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However, this introduces 3 more unknowns. For sensing 
coils that are not collinear or coplanar one can also 
write: 

cos a , cos a , + cos P , cos P , ♦ cos y , cos y , m k t (22 ) 

(cos p, cosy, - cosy, cos p,) 1 + (cosy, cosa, -cosa, cosy,) 1 + 

(cos a, cos P 2 - cos P, cosa,) J « *, (23) 

Equation (22) is derived from the definition from the 
vector dot product. Equation (23) is derived from the 
definition from the vector cross product. k 1 represents the 
cosine of the angle between the two sensor coil normals and 
k 2 represents the sine of the angle between the two sensor 
coil normals. When the sensor coils are orthogonal k t is 
identically 0 and k 2 is identically 1. This can be achieved 
readily by winding the sensor coils on a mandrel as shown 
in figure 6. When sensing coils are either coplanar or 
collinear k 1 is identically 1 and kj is identically 0. It 
should be noted again that when the sensing coil normals 
are collinear the system is only a 5DOF measurement system. 

Since the 2 sensing coils are not at the same position 
(even if they were wound around one another they might not 
be perfectly collocated) , the following substitutions 
(analogous to those for displacing the generating coils) 
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are required to reference the sensor coils to a common 
sensing point: 



x,«x-r,cos (a, ) (24) 
y^y-^cos^,) (25) 
2 1 -z-r,cos(Y, ) (26) 
and 

x 2 =x-r 2 cos(a 2 ) (27) 
y 2 =y-r 2 cos (0 2 ) (28) 
z 2 =z-r 2 cos(y 2 ) (29) 

where r, and r 2 represent the distance between the center of 
each sensor coil and the point x, y, z (see figure 7) . 
Now we can write: 



r j - - ' 2 J • if - T^F* it - » 2 ; oi i 

There are now ten equations (m=10) to choose from to solve 
for the 9 unknowns (n=9) . of course, one of the coil 
centers can be taken as the common reference point. This 
would eliminate one equation (m=9) . 

The direction cosines for an orthogonal set of sensor 
coils can be represented in terms of azimuth, elevation and 
roll (R) . For orthogonally placed sensing coils equation 

- 32 - 



WO 99/32033 



PCT/US98/13431 



(20) and (21) be rewritten as: 

(cos A cos E) 2 + (sin A cos £)' + (- sin £ )* -1 = 0 (32) 

(- cosR sin A + siart sin£cos^)* + (cos/? cos/4 + s in R sinEsin^) 1 

+ (sin*cos£) 2 -l=0 (33) 

This formulation is now also sufficient to solve the 
problem. Likewise # equation (22) and (23) can be rewritten 
5 in terms of azimuth, elevation and roll. This would 

provide additional equations for a non-linear least squares 
solution. 

2 Generating Coils, 3 Sensing Coils 

2 generating coils (p=2) can be placed at arbitrary 
positions and orientations as long as they are not 
collocated and oriented the same way. The coils should be 
arranged such that symmetries are minimized to prevent 
places having the same B-fields for different positions and 
orientations. The 2 generating coils may be coplanar or 
collinear or be a set of collocated orthogonal coils or 
anything else. Such different coil configurations are 
shown in FIG. 5. 

The 3 sensing coils (q=3) can be orthogonal, coplanar 
or collinear or some other combination as long as they are 
not collocated and oriented the same way. There are now 2 
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B-field equations per sensing coil for a total of 6 
equations. Either 5DOF or 6DOF systems can be constructed 
depending on whether the sensing coils are collinear, 
coplanar or whether they are oriented differently. The 
5 unknowns are the x, y and z position of the sensing coils 
and the 9 direction cosines. Expressed in terms of equation 
(11) there are 12 unknowns. 

Six or more additional equations are needed to solve 
for position and orientation. Three obvious equations are: 

io cos a a,+cos J P, +cos , y,-l»o o«) 

cos 2 a 2 +cos J p 2 +cos J Yj-l =0 ( 35 ) 

cos 2 a 3 + cos 2 p 3 +cos 1 y, -1 = 0 (36) 

However, this introduces 6 more unknowns. For sensing 
coils that are not collinear or coplanar one can also 
write: 



cos a, cos Oj + cos p, cos P, + cos y , cos y, - *, (37) 
cos a, cos a, + cos 0, cos p, + cos y , cos y, = *, ( 38) 
cos a , cos a , + cos p , cos p , + cos y , cos y , = A, (39) 

(cos p , cos y , - cos y , cos p, ) J + (cos y , cos a, - cos a , cos y , )* +■ 
(cosa.cospj-cosp.cosa,)' = k, (40) 
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(cos p | cos y 3 — cos T , cos P ,) 2 + (cos y , cos a , - cos a, cos y ,) a + 
(cos a, cos p, - cos p, cosa,) 2 = *, 



(cos p 2 cosy, -cosy, cos P,) 1 + (cos y x cosa, - cos a, cos y,) 1 + 
(cos a 2 cos p, - cos p 2 cosa,) 2 = k< (42 ) 



Equations (37-39) are derived from the definition from the 
vector dot product. Equations (40-42) are derived from the 
definition from the vector cross product. kj-k 3 represents 
the cosine of the angle between two sensor coil normals and 
k 4 -k 6 represents the sine of the angle between two sensor 
coil normals. When the sensor coils are orthogonal k t -k 3 
are identically 0 and k 4 -k 6 are identically 1. This can be 
achieved readily by winding the sensor coils on a mandrel 
as shown in figure 6. When sensing coils are either all 
coplanar or col linear k,-k 3 are identically 1 and k 4 -k 6 are 
identically <)• It should be noted again that when the 
sensing coils are all collinear the system is only a 5D0F 
measurement system. Other combinations such as 2 sensing 
coils coplanar and orthogonal to the third are also 
possible. 

Since the 3 sensing coils are not at the same position 
(even if they were wound around one another they might not 
be perfectly collocated) , the following substitutions 
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(analogous to those for displacing the generating coils) 
are required to reference the sensor coils to a common 
sensing point: 



10 



x 1 =x-r 1 cos(o, 


1 (43) 


yj-y-rjcostf, J 


(44) 


z i =z - r i c °s(Y, ) 


(45) 


and 




x 2 =x-r 2 cos (a 2 ) 


(46) 


y 2 =y-r 2 cos(/9 2 ) 


(47) 


2 2 = 2 -r 2 cos(y 2 ) 


(48) 


and 




x 3 =x-r 3 cos (Oj ) 


(49) 


yj=y-r 3 cos(0 3 ) 


(50) 


z 3 =z-r 3 cos(yj ) 


(51) 



15 



where r v r 2 and r 3 represent the distance between the center 
of each sensor coil and the point x, y, z. 
Now we can write: 



(52) 



(53) 



(54) 
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There are now 18 equations (m=18) to choose from to solve 
for the 12 unknowns (n=12) . Of course, one of the coll 
centers can be taken as the common reference point. This 
would eliminate one equation (m=17) . 

The direction cosines for an orthogonal set of sensing 
coils can be represented by azimuth, elevation and roll. 
For orthogonally placed sensing coils equations (31-33) can 
be replaced by 

(cos A cos E) 1 + (sin A cos E) 1 + (- sin E? - i = 0 (55) 

(-cos^siii/4 +sin^sin£cosy4)* +(cQ$RcosA+sinRsinEsinAf 
4-(sin/*cos£j l -l«0 (56) 

(sini?sin/4+cos^sin£ cos/4) 7 +(-sin/?cos/4+cos/?sin£sin^) 1 
+(cos/?cos£) 1 -1 = 0 (57) 



This formulation is also sufficient to solve the problem. 
Likewise, equations (37-42) can be derived in terms of 
azimuth, elevation and roll. This would provide additional 
equations for a non-linear least squares solution. 
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Obviously one can keep this process of generating 
configurations up indefinitely. 



Quality of Solution 

Ideally, the solution obtained fro* the chosen P & o process 
should, when inserted back into the defining equations 
(i.e., equation 11), yield zero. This is not likely to 
happen because of: 

Dthe finite resolution of the numerical representation 
within the computer, 

2) the finite resolution of the measured signals, 

3) environmental noise, 

4) imprecision of the system constants, 

5) distortion of the magnetic fields due to nearby 
conducting or ferromagnetic materials, etc. 

How accurately the P*o solution fits the equations is 
measured by the norm of the residuals. The norm of the 
residuals is defined as the square root of the sum of the 
squares of the values of the equations (11) evaluated at 
the P&o solution. The norm is what Hansen (4622644) calls 
the noise error vector. 

in practice a characterized system will produce a 
small but finite norm in an environment that contains no 
magnetic distorters (highly conductive materials and or 
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ferromagnetic materials) after the P&O process has 
converged. This value can be determined experimentally 
during system characterization in a distorter free 
environment and would, with careful design, cover the first 
4 items listed above. Therefore, major increases in the 
norm can be looked upon as a measure of the decrease in 
quality of the P&O solution due to distorters within a 
different environment, e.g., a catheter laboratory. This 
implies that the norm can be used to qualify whether the 
P&O solution is accurate enough for a particular 
application. By noting the range of acceptable norms in a 
particular environment, a threshold can be set by which the 
tracking quality can be monitored. When the norm is seen 
to exceed this threshold due to unexpected environmental 
noise or the inadvertent introduction of distorting 
materials a warning can be given such that an inaccurate 
P&O solution can be noted for evaluation. A flowchart for 
this process is shown in figure 15. 

Another use of the norm is in the compensation of 
magnetic field distortion caused by conductive and 
ferromagnetic materials. Assume that the norm is 
sufficiently small in a magnetically clean (no distorters) 
environment. Now the system is started at either a known 
P&O in the vicinity of distorters or in a region that is 
free of distorters. If started at a known P&O, then the 
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residuals, making up the norm are a measure of the 
distortion. These values can be filtered to eliminate 
noise impulses and subtracted from the next cycles P&o 
determination. If starting in a region without distortion 
the residuals are already small. As distorters are now 
introduced the magnitude of the residuals will increase. 
The residuals can again be filtered to eliminate noise 
impulses and subtracted from the next cycles P&o 
determination. A flowchart of this process is shown in 
figure 16. 



Starting Point 
A number of ways of providing an initial starting 
point are possible. One method is to perform a global 
search of the n dimensional space subject to the physical 
constraints of the application. This can be a very 
structured search of the n dimensional space, made by 
subdividing the space into an n dimensional grid and 
visiting every point on the grid. The subset of grid 
points with an acceptably small norm would then be tested 
in the P&o process until the global solution that best 
minimizes the norm is found. Another method is to perform 
an n dimensional random search of the space. The random 
points are tested in the P&o process until the best 
solution (minimal norm) is found. There are also systematic 
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methods as described by Aird, T. and Rice, J. in 
■Systematic Search in High Dimensional Sets,* SIAM Journal 
Numerical Analysis, Vol. 14 , No* 2, April 1977. Again, the 
points are tested in the P&O process until the best 
solution (minimal norm) is found. All of these methods 
suffer from the large quantity of points that need to be 
evaluated and are therefore expensive in terms of time. 
Hansen (4622644) mentions this approach. 

Another starting method is to use global methods such 
as bisection or homotopy methods. While these methods are 
expensive in terms of time they can usually find the 
correct P&O solution. 

Another way is to place the sensing coil(s) at a known 
P&O, which is noted in Dumoulin et al. (5377678). The 
process is then started with this P&O as it' s initial 
guess. This guarantees almost immediate convergence. This 
method can be extended to cover a final gain calibration, 
e.g., to correct for variations in the effective area (X q ) 
of the catheter sensor. Since the P£0 of the sensor is 
known by virtue of its placement, and, assuming all other 
characterizations are known, the remaining error in the 
norm can be assigned to X q . This amounts to adjusting A q to 
minimize the norm at the known placement and then using 
that value until a new catheter sensor is used. A 
flowchart of this process is shown in figure 17. 
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Once the process has converged correctly, that P40 
solution becomes the starting point for the next P&o 
measurement cycle. As long as drastic changes in actual 
P40 don' t occur between sensing coil measurements this is 
a valid and simple method for obtaining the processes next 
starting point. This also implies that the fastest 
processes for determining P &0 and the fastest sensed signal 
collection should be used to minimize the time between 
successive P&0 computations. This minimizes the sampling 
time and therefore the amount of movement the sensing 
coil(s) experiences. That way the starting point is very 
close to the correct P&o solution. 

Starting Point for 1 Sensing Coil 
As mentioned above, there are a number of ways of 
determining an initial starting point for the process. 
Under certain geometries and configurations the use of the 
infinitesimal dipole model can yield an approximate 
starting point. The use of the infinitesimal dipole model 
for determining a starting point does not invalidate the 
results disclosed in this patent. It is just another tool 
for determining a starting point. The following 
discussions are for circular coils. other numerical 
approximations could be determined by those skilled in the 
art for other geometries and configurations. Specifically, 
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a method for determining a starting point when using a 
single sensing coil can be realized if the generating coils 
are arranged as 2 (or more) sets of approximately 
collocated orthogonal triaxial coils or 3 (or more) sets of 
5 collocated orthogonal biaxial coils. This can be derived 
by utilizing the approximate B-field equations for coils in 
the y-z, x-z and x-y plane, respectively. 



h W^'.(( 2 x» -/-!•>* + 3<y3, *3«S ( ) (58) 

S _ l»± ln 2Lbxya t +(2/ -x % -x'fr, +3^5.) (59) 
' 4nr* 

S> ^ m b^SL(, aSt + »*3,+(2z> (60) 



where /x 0 is the permeability of free space; a is the coil 
10 radius; I is the current flowing through the coil; a K , a y 
and a z represent the unit vectors in Cartesian coordinates; 
r is the distance between the generating and sensing coil; 
and k represents the number of turns comprising the coil. 
Starting Point with 2 Triaxial Generating Coil Sets 
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Assume A p , A q , I p and a p are the same for all coils 
within a triaxial set and that the sets are collocated and 
orthogonal. Now assume that the sensing coil is oriented 
such that cos(o)=i, cos(/»)=o and cos(y)=o. Based on 
equation (11) the induced voltages for the triaxial coil 
set at the origin are: 



v ^'W^'V -/'-«') (6D 
v, ., - ' ' ' (3iy) (62) 

»u - r-i 



4«r* 



(63) 



Squaring and summing and taking the square root of the 
equations one gets: 

Using the binomial expansion on the r^x 2 term (assuming 
r 2 >3x 2 ) and simplifying equation (64) you obtain 
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1.1 



r j.i 



4nr» 



(65) 



This can be solved for r (now referred to as r 0 ) . Likewise, 
a r y can be determined for the other triad of generating 
coils. Referring to figure 8 one can see that knowing the 
r from both sets of triads defines a circle in the x-z 
plane with radius r e and centered at (0, y 0 , 0). This locus 
of points defines a fixed y coordinate and x and z 
coordinate pairs that lie on the circle 

rj » x 1 + 2 s (66) 



10 Points on the circle can now be tested in the P&O process. 

Starting Point with 3 Biaxial Generating Coil Sets 

Assume k pt X q , I p and a p are the same for all coils 
within a biaxial generating set and that they are 
collocated and orthogonal. Place coils in the y-z and x-z 
15 planes at the origin. Place this same configuration on the 
y-axis. Place a biaxial set on the x-axis, but this time 
place the coils in the x-y and x-z plane (see figure 9) . 
Now assume that the 1 sensing coil is oriented such that 
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cos(a)=o, cos(/J)=o and cos( Y )=l. Based on equation (11) 
the induced voltages due to the set at the origin are: 

The ratio of these two equations defines a line in the 
x-y plane 

P = f (69) 

Likewise, another line in the x-y plane can be derived for 
the biaxial set located on the y-axis. The intersection of 
these two lines yields a solution point x,y. This can be 
substituted into 

1 11 (70) 



4nr» Kiyt) 



which is the ratio of the voltages for the coils on the x- 
axis. The x and y values are substituted into this 
equation to solve for z. 
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Starting Point for 2 Sensing Coils 

The general comments from the previous section apply 
here also. A specific example for the 2 generating coils 
and 2 sensing coils starting point will be derived in a 
similar manner to that discussed in Starting Point with 2 
Triaxial Generating Coils, above. 

Starting Point with 1 Biaxial Generating Coil Set 

Assume JL, A„, I„ and a ft are the same for all coils 

p' q' P P 

within the biaxial generating set and that it is collocated 
and orthogonal. Place the coils in the y-z and x-z planes 
at the origin. Now assume that one sensing coil is 
oriented such that cos(a)=0, cos(/J)=0 and cos(y)=1 and the 
other sensing coil is oriented such that cos(a)=0, cos(/9)=0 
and cos(y)=0 (see figure 10). Based on equation (11) the 
induced voltages at the first sensing coil due to the 
generating coils at the origin are: 

*u a It i ( 3 *0 (71 > 
4nr 

= ' ' * \ 1 (3>*) <"> 
4nr 

The ratio of these two equations defines a line in the 
x-y plane. 
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i.i 



"7 <73> 



Likewise, the induced voltages at the second sensing coil 
due to the generating coils are: 



^.i.- 4nr» V* X y ) 



'j.i 



(74) 



5 Substituting equation (73) into (74) yields a curve, the 
points of which can now be tested in the P&O process. 

Starting Point with 1 Triaxial Generating Coil Set 

Assume k^, A„, r and a ft are the same for all coils 
within the triaxial generating set and that it is 
10 collocated and orthogonal. Place the coils in the x-y, y-z 
and x-z planes at the origin. Now assume that one sensing 
coil is oriented such the cos(o)=l, cos (6) and cos(y)=0 and 
the other sensing coil is oriented such that cos(a)=0, 
cos(6)=l and cosfYj^O. Based on equation (11) the induced 
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voltages at the first sensing coil due to the generating 
coils at the origin are: 



v,,, - ' ' fry) <w 



w * - (77) 



Voltages at the second sensing coil due to the generating 
coils at the origin are: 



(78) 



Squaring and summing and taking the square root of the 
equations one gets: 

Using the binomial expansion and simplifying equation 
(81) you obtain 
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4< 



V^, + V* , + v.* , + V,\ + V' , + V, 



"a.i 



'j.i 



1 J2 



4nr* 



(82) 



This can be solved for r. Equation (82) and ratios of v, ,/ 

v 3 , and v t 2 / v 3 2 can now be solved for x, y and z. 

Similar equations can be derived for other multiple 

sensor coil starting points. 

Preferred Embodiments 
The equations given above provide numerous different 
ways to set up a P&O system. It will be considered in the 
preferred embodiments that transducers (e.g., coils) used 
for generating information (e.g., electromagnetic fields) 
are external to the patient and the transducers used for 
sensing are in the catheter. However, one could easily 
reverse that so that generating transducers are in the 
catheter and sensing transducers are external to the 
patient and in the locator. Note that the catheter could 
be more broadly considered a remote device for which one 
wants to determine location parameters. As used herein, 
location parameters shall include position and orientation 
information relative to the catheter (or other remote 
device) . 
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The discussion will use a first example with eight 
coils used for field generation and one coil used for 
sensing. A second example with three generating coils and 
two sensing coils will then be discussed. These are the 
5 preferred embodiments when high accuracy and large regions 
of convergence are required. A third example will use four 
coils for field generation and one coil for sensing. A 
fourth example will use two coils for field generation and 
two coils for sensing. These last two examples are the 

10 preferred embodiments when low accuracy and limited regions 

of convergence are required. As mentioned previously, an 
important feature of the invention is the ability to 
provide 5 or 6 degrees of freedom location parameters when 
p*q=4. However, extra transducers may also be used to 

15 enhance the accuracy of the solution and the region of 
convergence . 

5DOF System 

The functional diagram shown in Figure 11 illustrates 
the preferred high accuracy 5DOF embodiment. A source of 
20 electromagnetic fields is generally illustrated at 10. The 
source of electromagnetic fields includes a plurality of 
field-generation means such as generator antennas or coils 
11 through 18 • While it is known that the minimum 
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configuration requires only 4 B-field generators, it has 
been found that 8 B-field generators provide a greater 
region. of convergence and a more accurate P40 solution. B- 
field generators placed as shown can utilize the fastest 
process by eliminating the possibility of rank deficiency 
in the P&o process. it is not necessary that there be 
precisely two B-field generators in a set or that there be 
eight B-field generators, in the present embodiment, the 
system operates in the near field. Amplifier 31 18 
establish quasi-stationary magnetic fields. Quasi- 
stationary magnetic fields are low frequency fields that 
vary so slowly with time that radiation effects are 
negligible. Although the present embodiment operates in the 
near field with coils, it should be appreciated that other 
embodiments may operate in the far field with other 
transmitting means, and rotating magnets can be substituted 
for the B-field generator coils. 

A signal generator is provided for applying electrical 
signals to B-field generator coils 11 through 18 to 
generate a plurality of low frequency electromagnetic 
field. Any frequency from zero to several hundred kHz may 
be used. The signals are multiplexed so that the fields 
generated by each of the coils are distinguishable. The 
signal generator is comprised of a digital to analog (D/A) 
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converter 41, a multiplexer 42 and driving amplifiers 71-78 
for supplying power to each of the B-field generator coils 
11 through 18 • The D/A 41 is fed a digitized waveform, 
typically sinusoidal, from the processor 50, which performs 
5 the P&O calculation process. It will be appreciated that 
the digitized waveform can be adjusted in gain by processor 
50 to affect controlled amplitude variations in the 
generated B-field. This would typically be used to 
compensate for large variations in sensed signal due to a 

10 large range of motion of the sensor, thus forming part of 
an automatic gain control (AGC) stage controlled by 
processor 50. The multiplexer comprises a conventional 
circuit for time division multiplexing the signals applied 
to each of the B-field generating coils. However, it will 

15 be appreciated by those skilled in the art that any one of 
a number of suitable multiplexing techniques may be 
employed, including time division multiplexing, frequency 
division multiplexing, or phase multiplexing and a suitable 
driver for implementing such a multiplexing technique will 

20 be provided. 

A sensor 20 receives the electromagnetic fields 
generated by the source 10. The sensor is a coil but other 
technologies such as flux gate, Hall effect, and magneto- 
resistive devices may be used. The sensing device must 
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produce an output proportional to the magnitude of the 
magnetic field and the cosine of the included angle between 
the directions of the magnetic field and the sensing device 
axis. For a coil the direction of the coil is 
perpendicular to the plane of the coil with the well-known 
right-hand-rule determining the sense. The output of the 
sensor coil is input into an amplifier 31. signal processor 
32 then processes the amplified signal. The purpose of the 
signal processor 32 is to reduce out of band signals from 
reaching the process processing stage (where P*o is 
calculated) and affecting the accuracy of the P&0 solution, 
in the preferred embodiment it is a simple low pass filter 
comprised of a resistor and a capacitor. It should be 
appreciated that any form of signal processing including 
none, passive, active analog and switched capacitor could 
be used depending on the application, m addition, signal 
processor 32 could contain adjustable gain stages to affect 
controlled amplification of the sensed sensor signal. This 
would form part of an AGC stage controlled by processor 50. 
An analog to digital conversion is then provided at 33. 

It should also be appreciated that any number of 
sensors 20 could be tracked using generating coils io, 
drivers 43 through 48, multiplexer 42, D/A 41 and processor 
50. Additional sensors could be multiplexed into amp 31 
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and then.be processed by signal processor 32, A/D 33 and 
processor 50. Alternately, multiple input channels 
consisting of sensor 20, amplifier 31, signal processor 32, 
A/D 33 could be multiplexed into processor 50. Those 
skilled in the art will obviously note other architectures. 

The signal from A/D 3 3 is then input into a suitable 
processor 50. The processor may be one or more actual 
processors optimized for running various parts of the 
system. A coherent detection is then performed within the 
processor between the digitized signal and a suitably phase 
shifted version of the digitized waveforms that drove D/A 
41. It should be appreciated that other signal extraction 
methods could be performed in place of or in addition to 
coherent detection, including FFTs, auto and cross- 
correlation, FIR and IIR filtering, etc. These signal 
processing methods can also be implemented outside the 
processor and within signal processor 32. Processor 50 
now determines the P&O of sensor 20. Processor 50 controls 
signal generation and detection, multiplexing and I/O with 
other systems. 
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As shown, the sensing coil 20 is disposed in the 
distal end of a catheter 60 having an optical fiber 62 
extending lengthwise therein. The axial direction of the 
cylindrical sensing coil 20 is preferably collinear with 
(or at least parallel to) the axis of the catheter 60. The 
portion of the catheter 60 housing the coil 20 is 
sufficiently rigid that coil 20 maintains its cylindrical 
configuration. Except for the details of the location 
determination arrangement described herein, the catheter 
may be otherwise constructed as described in the pending 
U.S. Patent application serial no. « fill this in» filed 
«fill in date» in the name of Abela and Bowden. That 
application is assigned to the assignee of the present 
application and is hereby incorporated by reference. 

The functional diagram shown in Figure 12 illustrates 
the preferred 5DOF embodiment when lower accuracy is 
acceptable. it is similar to the preferred embodiment 
illustrated in FIG. 11 with the exception that only 4 
generating coils u, 12, 13, and 14 and four drivers 71, 
72, 73, and 74 are required. Therefore, a detailed 
description of its operation and construction is not 
required. 
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In the discussion that follows, sensing coils of FIGS. 
13 and 14 may also be located in a catheter, although the 
catheter is not shown for ease of illustration. 

6DOF System 



5 The functional diagram shown in Figure 13 illustrates 

the preferred high accuracy 6DOF embodiment. The 
components of this embodiment have the same last two digits 
as the corresponding component of the FIG. 11 embodiment. 
A source of electromagnetic fields is generally illustrated 

10 at 110. The source of electromagnetic fields includes a 
plurality of field-generation means such as generator 
antennas 111, 112 and 113. It is not necessary that there 
be precisely three B-field generators in a set or that 
there be one B-field generator set. In the present 

15 embodiment, the system operates in the near field. 

Antennas 111 through 113 establish quasi-stationary 
magnetic fields. Quasi-stationary magnetic fields are low 
frequency fields that vary so slowly with time that 
radiation effects are negligible. Although the present 

20 embodiment operates in the near field with coils, it should 
be appreciated that other embodiments may operate in the 
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*« field „ ith other transmitting ^ ^ 
»ag„ets can be substituted for the B-field generators. 

A sign.! generator is provided for applying electrical 
signals to B-field generator coils m through „, tg 
generate e plurality of low frequency electromagnetic 
field. Any f requ ency fro, 2 ero to several hundred m say 
be used. The signals are multiplexed so that the fields 
generated by each of the coils are distinguishable. The 
signal generator is comprised of a digital to analog (D/A) 
converter M1 . . multiplexer 142 and driving amplifiers 
143-145 for supplying power to each of the B-field 
generating coils m through 113. Tne D/A 141 ls f- , 
digitised waveform, typically sinusoidal. fro. the 
processor 150. It will be appreciated that the digitized 
waveform can be adjusted in gain by processor 50 to affect 
controlled amplitude variations in the generated B-field 
This would typically be used to compensate for large 
variations in sensed signal due to a large range o, motion 
f the sensor, thus forming part of an AGC stage controlled 
by processor 50. The multiplexer comprises a conventional 
~rcn.it for time division multiplexing the signals applied 
to each of the B-fie!d generating coils. However, it will 
be appreciated by those sxilled in the art that any one of 
e number of suitable multiplexing techniques may be 
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employed, including time division multiplexing, frequency 
division multiplexing, or phase multiplexing and a suitable 
driver for implementing such a multiplexing technique will 
be provided. 

5 A sensor 120 receives the electromagnetic fields 

generated by the source 110. The sensor is comprised of 
coils but other technologies such as flux gate, Hall 
effect, and magneto-resistive devices may be used. The 
components comprising the sensor must produce an output 

10 proportional to the magnitude of the magnetic field and the 
cosine of the included angle between the directions of the 
magnetic field and the component axis. The coil outputs 
are input into multiplexer 134 where each coil can be 
sampled in a time-multiplexed fashion. The selected coil 

15 signal is input into an amplifier 131. Signal processor 132 
then processes the amplified signal. The purpose of the 
signal processor 132 is to reduce out of band signals from 
reaching the processing stage (where P&O is calculated) and 
affecting the accuracy of the P&O solution. In the 

20 preferred embodiment it is a simple low pass filter 
comprised of a resistor and a capacitor. It should be 
appreciated that any form of signal processing including 
none, passive, active analog and switched capacitor could 
be used depending on the application. In addition, signal 
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processor 32 could contain adjustable gain stages to affect 
controlled amplification of the sensed sensor signal. This 
would ?orm part of an AGC stage controlled by processor 50. 
An analog to digital conversion is then provided at 133. 

It should also be appreciated that any number of 
sensors 120 could be tracked using generating coils no, 
drivers 143-145, multiplexer 142, D/A 141 and processor 
150. Additional sensors could multiplexed into amp 131 via 
multiplexer 134 and then be processed by signal processor 
132, A/D 33 and processor 150. or multiple input channels 
consisting of additional sensors 120, multiplexer 134, amp 
131, signal processor 132, A/D 133 could be multiplexed 
into processor 150. Those skilled in the art will 
obviously note other architectures. 

The signal from A/D 133 is then input into a suitable 
processor 150. The processor may be one or more actual 
processors optimized for running various parts of the 
system, a coherent detection is then performed within the 
processor between the digitized signal and a suitably phase 
shifted version of the digitized waveforms that drove D/A 
141. It should be appreciated that other signal extraction 
methods could be performed in place of or in addition to 
coherent detection, including FFTs, auto and cross- 
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correlation, FIR and IIR filtering, etc. These signal 
processing methods can also be implemented outside the 
processor within signal processor 132. Processor 150 now 
determines the P&O of sensor 120. Processor 150 controls 
5 signal generation and detection, multiplexing and I/O with 
other systems. 

The functional diagram shown in Figure 14 illustrates 
the preferred 6DOF embodiment when lower accuracy is 
acceptable. It is similar to the preferred embodiment 
10 illustrated in FIG. 13 with the exception that only 2 
generating coils are required. The components of this 
embodiment are given the same numbers as the corresponding 
components of FIG. 13. 

System Parameter Determination 

The equations presented so far have been in terms of 
known parameters X p , A q , I p , a p , k, r (the distance between 
the center of each sensor coil and the point x, y, z in 
equations 24-31 and 43-54), and rotation matrices. These 
parameters need to be determined accurately before any P&O 
process is implemented. I p can be determined within the 
tracking system via specific measurement circuitry. The 
other parameters can be determined by running the P&O 



15 



20 
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process in reverse. 



specifics, if . sensor wlt „ ^ parMeters 
(obtained by other means) is placed at a altitude of Known 
P«0s, then the other parameters can be solved for using a 
least scares approach. The starting point is based on 
mechanical measurements of the B-,ield generator coil 
positions and orientations and manufacturing specifications 
of both generating and sensing coils. The data is most 
easily collected with a mechanic X -y- 2 translation table 
that moves a mechanized gimbal, such as schematically shown 
in PIG. 15. The gimbal may have the sensing coils on it 
and the generating coils fixed nearby or vice versa. The 
gimbal is designed to allow the sensing device to be placed 
at various orientations, o, course, the table and gimbal 
must be more accurate than the parameters you are trying to 
determine. It is assumed that the P40 of the sensing 
device within the gi^al is always Known. This alWs Me 
to calibrate the relation between the sensing 3n3 
generating devices such that the measurements made later on 
actual patients will be highly accurate. 

FIG. 16 shows a simplified block diagram of the 
quality of solution process used by the processors such as 
processor 50 of FIG. 11. Initially , . tracking ^ 
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would be evaluated by measuring the norms within a 
specified volume, as noted in box 200. This would be 
performed as part of the manufacturing process in a 
distortion free environment. The norm measures the residual 
5 error remaining after a proposed P&O solution is plugged 
back into the defining equations and indicates how accurate 
the P&O solution is. The worst case norm would be noted 
for the system and shipped with the system. It would also 
be used to set a related threshold level as shown in box 

10 202. Once the tracker is placed in use by a customer, all 
subsequent P&O solution norms calculated in box 204 would 
be compared to the threshold level. If the P&O solution 
norm exceeds the threshold level an indication would be 
made to the user signifying that the P&O solution is 

15 suspect, as shown in box 206. This process continues 
indefinitely, comparing each new solution norm to the 
threshold. This process thereby alerts the user to 
problems within the environment. Once the environmental 
problem is corrected the indicator is automatically reset 

20 by the tracking system. Alternately, the norms could be 
monitored over the volume of interest once the tracking 
system was placed in operation. This process would 
typically be done once the system was installed in its 
final tracking environment. This way the worst case norm 

25 represents the specifics of the users environment. In 
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either case the introduction of distorters to the tracking 
environment would increase the norm above the threshold 
level, indicating that the solution is "not as good as it 



once was. 



FIG. 17 shows a simplified block diagram of a process 
used with the present invention in the compensation of 
magnetic field distortion caused by conductive materials. 
Assume that the norm is sufficiently small in a 
magnetically clean (no distortion) environment. Now the 
system is started at either a known P&o in the vicinity of 
distorters or in a region that is free of distorters as 
shown in block 210. if started at a known P&o, then the 
residuals making up the norm are a measure of the 
distortion. These values can be low pass filtered to 
eliminate noise impulses and provide short-term memory as 
shown in block 212. If starting in a region without 
distortion the residuals are already small. As distorters 
are now introduced the magnitude of the residuals will 
increase. The residuals can again be low pass filtered to 
eliminate noise impulses as noted in block 212. In either 
case the filtered residuals are subtracted from the 
defining equations as shown in block 214. The next 
computation cycle of P&o determination then uses the 
adjusted equations as noted in block 216. This process 
25 continues indefinitely. 
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FIG. 18 shows a sensor calibration process that may be 
performed using any of the tracking systems of FIG. 11-14 • 
The assumption here is that the rest of the system has 
already been characterized using the apparatus shown in 
5 FIG. 15 and only the sensor needs calibration. An example 

of this would be when a new sensor is used with a 
calibrated tracking system. At block 218, the sensor is 
placed at a known position and orientation. This could be 
a special mount or reference point that is fixed in 

10 relation to the B-field generating coils. The P&O data is 
calculated at block 220 assuming a nominal sensor gain. At 
block 222 , deviations between the known P&O and calculated 
P&O are assigned to errors from the sensor gain . A new 
sensor gain is established in block 224. The new sensor 

15 gain may then be used for actual patient measurements. 

Disposable Sensor 

Turning now to FIGS. 19-23, a specific type of 
disposable sensing coil that may be used with the present 
invention is discussed. The sensing coil 300 is disposable 
2 0 as that term is used in the medical field where certain 
medical items are designed with disposal in mind. These 
include items that cannot easily be sterilized and items 
that are not cost-effectively sterilized. 
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The 5DOF tracking implementation provides the ability 
to construct inexpensive sensing coils for tracking 
exterior portions of the body or exposed interior portions. 
A single sensing coil 300 can be built by laying out a 
flexible circuit board 302 composed only of a circular 
track pattern 304 (PIG. 19). other patterns are possible 
including a spiral pattern 306 (FIG. 2) . For all of the 
sensors discussed with reference to FIGS. 19-23, an 
adhesive patch backing 308 (shown in FIG. 19 only) is used 
to adhere the sensor to part of a patient. Leads would be 
integrally attached and terminated with a simple 2 contact 
connector (FIG. 19) . 

A multiple layered flexible circuit board 310 allows 
the coil pattern to be extended in the z direction (FIG. 
21). Both of the aforementioned methods improve the coil 
sensitivity. The sensor board 310 has two boards 312 
adhered together (shown offset for ease of illustration) 
and two copper trace coils 314. 

Another possibility is to use flexible circuit board 
material 320 for the leads too (FIG. 22). a further 
possibility would be to incorporate a concentric coaxial 
snap connector 322 onto the flexible circuit for quick 
connect/disconnect capabilities (FIG. 23). 
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The flexible circuit board is mounted to an adhesive 
patch (PIG. 19). Both items can be produced in a sterile 
environment and transported in a suitable sterile package. 

Disposable Sensor Uses 

The disposable sensor can be used in many of the 
same applications as a non-disposable and various sizes 
could be constructed for different applications. Gait 
analysis setup would be facilitated because you would 
just place the patches at the point you want to monitor 
and you would be done. Velcro straps, large elastic 
bandages and other cumbersome attachment methods are 
presently used. The alignment of joints in arthroplasty 
procedures could be drastically improved by measuring the 
position and orientation of bone alignment before and 
during the procedure. In this case a bone nail or pin 
could be pushed through the center of the disposable 
sensor and attached at a suitable point of reference. 

Uses 

The remote object position and orientation 
determining system of the present invention has a wide 
variety of applications. For example, sensor 20 can be 
placed at the end of a catheter or flexible endoscope for 
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tracking the path of travel or the location of the 
catheter tip. This is especially useful in cardiac 
applications, e.g., for the treatment of cardiac 
arrhythmias, drug delivery and biopsy needle placement. 
Another use for sensor 20 is to embed the sensor into a 
contact lens for eye tracking. The eye tracking could 
also facilitate laser eye surgery by adjusting the aim of 
the laser during eye movements. Another application 
would place numerous sensor 20s around the fingers to 
provide a measure of finger movement for evaluating 
rehabilitation progress or as part of a virtual reality 
data glove, a device that allow hand movement to be 
translated into a virtual environment. Another 
application of sensor 20 would be to align prostheses 
during arthroplasty procedures. Sensor 20 can also be 
associated with the stylus of a three dimensional 
digitizer which is used to trace a physical model or the 
like and generate a digital database. The resulting 
database can then be used to generate a wide variety of 
computer generated images of the physical model, it 
should also be obvious that a 2 dimensional digitizer can 
be constructed using this technology, still another 
application would be as a stylus input for a Personal 
Digital Assistant (PDA) . 
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Sensor 120 can be used for forming three dimensional 
echocardiograms by tracking a handheld ultrasound scanner 
head. In still another application, sensor 120 can be 
associated with a particular body part for the purpose of 
conducting biomechanical studies. Still another 
application involves the monitoring of the body movements 
of an invalid for the purpose of creating a nonverbal 
communication system or providing a technique for 
remotely controlling various devices with nonverbal 
communicative body motion. Still another application is 
the production of data gloves , devices that allow hand 
movement to be translated into a virtual environment. 
Another application is the tracking of head movements for 
use in virtual reality headsets. Still another 
application is a 3 dimensional mouse for navigating in 
multidimensional databases or as a tool for use with 
VRML, a 3 dimensional file format popular on the 
Internet. 

Following is a list of documents providing 
mathematical and related background discussion and as 
referred to above by number: 
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art. In view of possible modifications, it will be 
appreciated that the scope of the present invention 



WO 99/32033 

PCT/US98/13431 



should be determined by reference to the claims appended 
hereto. 



- 74 - 



WO 99/32033 



PCT/US98/13431 



What is claimed is: 



1 1. A medical system for use on a patient 

2 comprising: 

3 a plurality of transducers operable to determine the 

4 location of a medical device in a patient, the 

5 transducers including one or more device transducers 

6 and one or more external transducers; 

7 a medical device for insertion into a patient, the 

8 medical device having the one or more device 

9 transducers ; 

10 an external locator operable outside and adjacent to 

11 the patient, the external locator having the one or 

12 more external transducers; 

13 a driver connected to sequentially energize one or 

14 more of the device transducers and/or the external 

15 transducers; 

16 a processor operably connected to receive sensor 

17 signals from one or more of the transducers, each 

18 sensor signal on a given one of the transducers 

19 depending on the energizing of another of the 
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20 

21 

22 

23 

24 

25 



26 validity. 



1 

2 



transducers, the processor operable to iteratively 
determining at least two location parameters of the 
medical device in the patient using a plurality of 
equations, which equations do not require a minimum 
separation distance from the energized transducer to 
the transducer generating the sensor signal for 



1 2. The medical system of claim 1 wherein the medical 

2 device is a catheter. 



3. The medical system of Claim 2 wherein the 

2 catheter is a laser catheter operable for endomyocardial 

3 revascularization. 



4. The medical system of claim 3 wherein each of 
the transducers is operable to generate and/or sense an 



3 electromagnetic field. 
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5. The medical system of Claim 4 wherein the 
processor determines n location parameters of the laser 
catheter and the number of transducers is less than n. 



6. The medical system of Claim 4 wherein there is 
only one device transducer and there are a plurality of 
external transducers. 



7. The medical system of Claim 4 wherein the 
external transducers generate electromagnetic fields and 
each device transducer senses the electromagnetic fields. 



8. The medical system of Claim 2 wherein each of 
the transducers is operable to generate and/or sense an 
electromagnetic field. 
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9. The medical system of Claim 8 wherein the 
processor determines n location parameters of the laser 
catheter and the number of transducers is less than n. 



10. The medical system of Claim 1 wherein the 
processor determines n location parameters of the laser 
catheter and the number of transducers is less than n. 



11. A system for determining the location of a 
remote device comprising: 

a plurality of transducers operable to determine the 
location of a remote device in an operation zone, 
the transducers including one or more device 
transducers and one or more locator transducers; 

a remote device for placement in an operation zone, 
the remote device having the one or more device 
transducers ; 

a locator operable away from, but within an 
interaction distance from, the remote device, the 
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12 locator having the one or more locator transducers; 

13 a driver connected to sequentially energize one or 

14 more of the device transducers and/or the locator 

15 transducers; 

16 a processor operably connected to receive sensor 

17 signals from one or more of the transducers , each 

18 sensor signal on a given one of the transducers 

19 depending on the energizing of another of the 

20 transducers, the processor operable to iteratively 

21 determining at least two location parameters of the 

22 remote device in the patient using a plurality of 

23 equations, which equations do not require a minimum 

24 separation distance from the energized transducer to 

25 the transducer generating the sensor signal for 

26 validity; and 

27 wherein each of the transducers generates and/or senses 

28 an electromagnetic field. 



1 12. The system of Claim 11 wherein the system is a 

2 medical system, the remote device is a medical device and 

3 the operation zone is within a patient. 
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13. The system of Claim 12 wherein the medical 
device is a laser catheter operable for endomyocardial 
revascularization. 



14. The system of claim 13 wherein the processor 
determines n location parameters of the laser catheter 
and the number of transducers is less than n. 



15. The system of claim 14 wherein n is at least 

five. 



16. The system of Claim 12 wherein there is only one 
device transducer and there are a plurality of locator 
transducers. 



17. The system of Claim 11 wherein there is only one 
device transducer and there are a plurality of locator 
transducers. 



- 80 - 



WO 99/32033 



PCT/US98/13431 



1 18. The system of Claim 11 wherein the locator 

2 transducers generate electromagnetic fields and each 

3 device transducer senses the electromagnetic fields. 

1 19. The system of Claim 11 wherein the processor 

2 determines n location parameters of the remote device and 

3 the number of transducers is less than n. 

1 20. The system of Claim 19 wherein n is at least 

2 five. 

1 21. A method of determining the location of a remote 

2 device comprising the steps of: 

3 using a system having a plurality of transducers 

4 operable to determine the location of a remote 

5 device, the transducers including one or more device 

6 transducers and one or more locator transducers; 

7 placing the remote device in an operation zone, the 
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8 remote device having the one or more device 

9 transducers; 

10 placing a locator away from, but within an 

11 interaction distance from, the remote device, the 
locator having the one or more locator transducers; 



12 



13 
14 

15 
16 
17 



sequentially energizing one or more of the device 
transducers and/or the locator transducers; 

generating sensor signals from one or more of the 
transducers, each sensor signal on a given one of 
the transducers depending on the energizing of 

18 another of the transducers; and 

19 supplying the sensor signals to a processor, the 

20 processor determining at least two location 
parameters of the remote device using a plurality of 
equations, which equations do not require a minimum 
separation distance from the energized transducer to 
the transducer generating the sensor signal for 

25 validity; and 



21 
22 
23 
24 



26 wherein the system is a medical system used on 

27 patient, the remote device is a medical device and the 

28 operation zone is within a patient. 
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1 22. The method of Claim 21 wherein each of the 

2 transducers generates and/or senses an electromagnetic 

3 field. 

1 23. The method of Claim 22 wherein each of the 

2 transducers is a coil. 

1 24. The method of Claim 21 wherein the processor 

2 determines a number N of locator parameters and there are 

3 N-l transducers. 

1 25. The method of Claim 21 wherein the transducers 

2 are selected from the group consisting of: coils, Hall- 

3 effect sensors, magnetometers , magnetoresistive devices, 

4 magnetoinductive devices, and galvanomagnetic devices. 

1 26. The method of Claim 21 wherein the medical 

2 device is a catheter. 
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27. The method of claim 21 wherein the medical 
device is a laser catheter. 



28. The method of Claim 27 further comprising the 
steps of: 

positioning the laser catheter within the heart of a 
patient using the determined location parameters to 
move the laser catheter to a desired location; and 

lasing from the laser catheter so as to perform 
endomyocardial revascularization. 



29. A method of determining the location of a remote 
device comprising the steps of: 

using a system having a plurality of transducers 
operable to determine the location of a remote 
device, the transducers including one or more device 
transducers and one or more locator transducers; 

placing the remote device in an operation zone, the 
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remote device having the one or more device 
transducers ; 

placing a locator away from, but within an 
interaction distance from, the remote device, the 
locator having the one or more locator transducers; 

sequentially energizing one or more of the device 
transducers and/or the locator transducers; 

generating sensor signals from one or more of the 
transducers, each sensor signal on a given one of 
the transducers depending on the energizing of 
another of the transducers; and 

supplying the sensor signals to a processor, the 
processor determining at least two location 
parameters of the remote device using a plurality of 
equations, which equations do not require a minimum 
separation distance from the energized transducer to 
the transducer generating the sensor signal for 
validity; and 

wherein each of the transducers generates and/or senses 
an electromagnetic field. 
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30. The method of claim „ vhereln the ^ ie a 

» system used to treat a patient> ^ ^ 

» i» « -did device end the operation zone is within . 
4 patient. 



31. The method of claim 29 wherein *-h„ 

« wnerein the processor 

determines the at ie.st two location parameters of the 
remote device by checking to insure that residue! errors 
are less than a given threshold. 



32. The method of claim 29 wherein the processor 
determines the at least two location parameters of the 
remote device by talcing into account previously 
determined distortions. 



33. The method of claim 32 further including the 
steps of, prior to the placement o, the remote device in 
*n operation zone, placing the remote device in a Known 
position relative to the locator, determining the 
deviation between sensed position data and known position 
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6 data, and adjusting for distortions based on the 

7 deviations such that the at least two location parameters 

8 of the remote device by take into account distortions. 
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